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ABSTRACT 

We measured the bar pattern speed, Q p , of the SBO galaxy NGC 1023 using 
the Tremaine-Weinberg (1984) method with stellar-absorption slit spectroscopy. The 
morphology and kinematics of the H I gas outside NGC 1023 suggest it suffered a 
tidal interaction, sometime in the past, with one of its dwarf companions. At present, 
however, the optical disc is relaxed. If the disc had been stabilized by a massive dark 
matter halo and formed its bar in the interaction, then the bar would have to be slow. 
We found fl p = 5.0 ± 1.8 km s _1 arcsec _1 , so that the bar ends near its co-rotation 
radius. It is therefore rotating rapidly and must have a maximum disc. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: haloes — galaxies: 
individual: NGC 1023 — galaxies: interactions — galaxies: kinematics and dynamics 
— galaxies: photometry 


1 INTRODUCTION 

Strong bars are seen in optical images of roughly 30% of 
all high surface brightness (HSB) disk galaxies (Sellwood 
& Wilkinson 1993) and this fraction rises to 50%-75% in 
the near IR (Knapen 1999; Knapen et al. 2000; Eskridge 
et al. 2000). Understanding the structure and dynamics of 
barred (SB) galaxies is, therefore, an issue of some impor¬ 
tance. The principal dynamical quantity for SB galaxies is 
the rotation frequency/pattern speed of the bar, fl p . This 
is usually parametrized by the distance-independent ratio 
TZ = Di/aB, where as is the semi-major axis of the bar 
and Dl is the radius to the Lagrangian point, where the 
gravitational and centrifugal forces cancel out in the bar’s 
rest frame. (The Lagrangian radius is therefore the general¬ 
ization to strong bars of the corotation radius.) Contopoulos 
(1980) argued that a self-consistent bar is possible only if 
TZ > 1. A bar is termed fast when 1.0 < TZ £ 1.4, while, for 
a larger value of TZ, a bar is said to be slow. 

A variety of methods have been used to attempt mea- 
surent of bar pattern speeds (see, for example, the review of 

* Based on observations made with the Italian Telescopio 
Nazionale Galileo (AOT-3, 3-06-119) operated on the island of La 
Palma by the Centro Galileo Galilei of the Consorzio Nazionale 
per l’Astronomia e l’Astrofisica, and the Jacobus Kapteyn Tele¬ 
scope operated by the Isaac Newton group at La Palma island 
at the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 
t email: debattisOastro .unibas. ch The observational evidence, therefore, favors fast bars. 


Elmegreen 1996). Most measurements of TZ rely on hydrody- 
namical simulations. These usually try to match the gas flow 
in the region of the bar, particularly at the shocks, which 
works because the location of the shocks depends on TZ, mov¬ 
ing further ahead of the bar as TZ increases. A bar needs to 
be fast for the shocks to remain curved with their concave 
sides towards the bar major axis, as is usually observed (van 
Albada & Sanders 1982; Athanassoula 1992). Detailed sim¬ 
ulations of gas flows in individual galaxies also result in fast 
bars; examples include: NGC 1365 (TZ = 1.3, Lindblad et al. 
1996), NGC 1300 (TZ = 1.3, Lindblad & Kristen 1996), and 
NGC 4123 (TZ = 1.2, Weiner et al. 2001). Hydrodynamical 
simulations can also recover TZ by matching morphological 
features in H I; some examples are: NGC 7479 (TZ = 1.22, 
Laine 1996), NGC 1073 (TZ = 1 - 1.2, England et al. 1990), 
NGC 3992 (TZ = 1, Hunter et al. 1989), and NGC 5850 
(TZ = 1.35, Aguerri et al. 2001). 

A direct method for measuring using a tracer popu¬ 
lation which satisfies continuity, was developed by Tremaine 
& Weinberg (1984). Since gas is subject to phase changes, it 
is not well-suited for this application. Old stellar populations 
in the absence of significant obscuration, on the other hand, 
are ideal for the Tremaine-Weinberg (TW) method. This has 
permitted application of the method to a small number of 
early type SB galaxies: NGC 936 (TZ = 1.4 ±0.3, Kent 1987 
and Merrifield & Kuijken 1995), NGC 4596 (TZ = 1.151q' 23 , 
Gerssen et al. 1998) and NGC 7079 (TZ = 0.9 ±0.15, Debat- 
tista & Williams 2001). 
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The perturbation theory calculations of Weinberg (1985) 
predicted that a fast bar would be slowed down rapidly 
in the presence of a massive dark matter (DM) halo. Such 
slow-down has been seen in various simulations (Sellwood 
1980; Little & Carlberg 1991; Hernquist & Weinberg 1992; 
Athanassoula 1996). The fully self-consistent, high resolu¬ 
tion IV-body simulations of Debattista & Sellwood (1998) 
also confirmed this prediction; however they showed that, 
for a maximum disc (here taken to mean a disc which domi¬ 
nates the rotation curve throughout the inner few disc scale- 
lengths, cf. van Albada & Sancisi 1986), a fast bar can sur¬ 
vive for a large fraction of a Hubble time. Subsequently, 
Tremaine & Ostriker (1999) suggested that bars manage to 
remain fast not because discs are maximal but rather be¬ 
cause the inner parts of DM halos are flattened and rapidly 
rotating. However, Debattista & Sellwood (2000) showed 
that rapid slow-down persists even then unless the halo an¬ 
gular momentum is very large relative to that of the disc. 
Thus they concluded that SB galaxies must be maximal, 
and argued that the same must be true for all high surface 
brightness disc galaxies. 

This conclusion rests on a small number of pattern 
speed measurements; in view of the fact that maximum 
discs are in conflict with the predictions of cold dark mat¬ 
ter (CDM) cosmologies ( e.g. Navarro et al. 1997), enlarging 
the sample of measured pattern speeds is desireable. In this 
paper, we report observations of NGC 1023, for which we 
applied the TW method. 

The rest of this paper is organized as follows. The TW 
method is described briefly in §2j. Then, in : jd] we give an 
overview of the previously known properties of NGC 1023. 
The photometric observations, reduction and results, includ¬ 
ing as, are presented in §^, while (jo] presents the spec¬ 
troscopic observations and results. We derive the rotation 
curve, corrected for the asymmetric drift, from which we 
obtain Dl. With these results at hand, we then apply the 
TW method in §|fij. We present our conclusions in » 7| 


2 THE TREMAINE-WEINBERG METHOD 


The TW method is contained in the following simple equa¬ 
tion: 


Dp sin i 


IZ Hr)fZ s ( x ’ Y ) v los (x,Y) dx dY 
fl 1 h(Y) fZ, £(X, Y) X dX dY 


(1) 


where i is the inclination of the galaxy, Vi os is the line-of- 
sight (minus systemic) velocity, E is the surface brightness, 
h is an arbitrary weighting function and (X, Y) are coor¬ 
dinates on the sky, centered on the galaxy, along the disc’s 
apparent major and minor axes, respectively. Slit observa¬ 
tions are equivalent to h(Y) oc 5(T — Y s u t ). Application of 
Eqn. |I] presents three main difficulties: 

(i) Centering errors. As was already recognized by 
Tremaine & Weinberg (1984), small errors in identifying ei¬ 
ther the center or the systemic velocity of the galaxy can sig¬ 
nificantly affect the value of D p obtained. Noting that most 
SB galaxies are nearly point-symmetric about their centers, 
Tremaine & Weinberg (1984) suggested using a weighting 
function which is odd in Y to counter this problem. Merri- 
field & Kuijken (1995) used a different, and somewhat bet¬ 


ter, strategy. Noting that the integrals in Eqn. hlrepresent 
luminosity-weighted averages, they rewrote Eqn7ll| as: 

DpSin^A — X c ) = V — V sys , (2) 

where V and X are the luminosity-weighted averages relative 
to an arbitrary frame, in which the galaxy center is (X c , Y c ) 
and the systemic velocity is V sys . Plotting V against X then 
gives a straight line with slope DpSini. With this approach, 
the problem of centering errors becomes one of fixing an 
arbitrary reference position and velocity frame common to 
all the slits, which in general is much easier to achieve. 

(ii) Small signal-to-noise ratios. V and X measure dif¬ 
ferences across X = 0 and are, therefore, susceptible to 
noise. This is particularly true for V, for which the non- 
axisymmetric part of V) os is not much larger than typical 
measurement errors. Merrifield & Kuijken (1995) overcame 
this problem by projecting their slit spectra along the spa¬ 
tial direction thereby increasing the signal-to-noise ratio of 
V significantly. This projection amounts to carrying out the 
required velocity integral directly in photon space. The sig¬ 
nal in X is typically much better constrained, since bars 
tend to be bright. 

(iii) Sensitivity to errors in the disc’s position angle. 
This was first recognized by Debattista & Williams (2002). 
They used two-dimensional Fabry-Perot absorption-line 
spectroscopy of NGC 7079 to show that small errors in the 
derived position angle (PA) of the disc translate into large 
errors in D p . For example, for NGC 7079, they found that an 
error of as little as 5° in the disc PA would result in an error 
of 100% in D p . Since the uncertainties in published values of 
PA for disc galaxies are often at the 5° level, each applica¬ 
tion of the TW method requires careful measurement of the 
disc PA. Highly inclined galaxies and galaxies in which the 
bar is at about 45° to the disc major axis, are less sensitive 
to errors in disc PA. We note that past applications of the 
TW method using slits have been to well studied galaxies 
(Merrifield & Kuijken 1995; Gerssen et al. 1998), so that it is 
very unlikely that large errors of this type were introduced. 

In this work, we have followed the prescriptions of Mer¬ 
rifield & Kuijken (1995) to deal with the centering and 
signal-to-noise problems. We dealt with the sensitivity to 
errors in disc PA by obtaining deep surface photometry of 
NGC 1023 in advance of the spectroscopy. 


3 GLOBAL PROPERTIES OF NGC 1023 

NGC 1023 (UGC 2154) is a highly-inclined lenticular galaxy 
classified as SB0 by Nilson (1973), as SB0i(5) by Sandage 
& Tammann (1981), and as SB0 _ (rs) by de Vaucouleurs 
et al. (1991, hereafter RC3). Its total B-band magnitude 
is Bt = 10.35 (RC3), which, after correcting for inclina¬ 
tion and extinction, corresponds to Mg T = —19.96 for an 
adopted distance of 10.2 Mpc (Faber et al. 1997; Ho = 80 
km s -1 Mpc -1 ). 

NGC 1023 is the brightest member of the LGG 70 
group (Garcia 1993). Its closest companion, designated NGC 
1023A by Hart et al. (1980), is a low luminosity condensation 
located near the eastern end of the major axis. The centers of 
the two objects are separated by an angular distance of 2'.7 
(RC3), corresponding to a projected linear distance of about 
8 kpc at 10.2 Mpc. The proximity of NGC 1023A to NGC 
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1023 led Arp (1966) to include NGC 1023 in the section of 
his Atlas of Peculiar Galaxies (Plate 135) devoted to E-like 
galaxies with nearby fragments. NGC 1023A was recognized 
as an individual galaxy by Barbon & Capaccioli (1975) and 
classified as Magellanic irregular or late-type dwarf galaxy 
by Capaccioli, Lorenz & Afanasiev (1986). High resolution 
radio observations by Sancisi et al. (1984) found a complex 
structure and kinematics for the neutral hydrogen. Most of 
it was shown to be outside of the galaxy and consists of 
tails and bridges connecting NGC 1023 to three compan¬ 
ions, including NGC 1023A. This morphology suggests an 
interaction with NGC 1023A, a conclusion also reached by 
Capaccioli et al. (1986). Although it is rich in neutral hydro¬ 
gen, NGC 1023 is undetected in ionized (Pogge & Eskridge 
1993; Sanchez-Portal et al. 2000) and molecular gas (Sofue 
et al. 1993; Taniguchi et al. 1994), making application of the 
TW method to it possible. 

Surface photometry of NGC 1023 has been obtained in 
several optical (e.g., Barbon & Capaccioli 1975; Kormendy 
1985; Lauer et al. 1995; Sanchez-Portal et al. 2000; Bower 
et al. 2001) and near-infrared bands (Mollenhoff & Heidt 
2001). Detailed stellar kinematics have been measured by 
Simien & Prugniel (1997, PA = 87°, 177°), Neistein et al. 
(1999, PA = 87°), and recently by Bower et al. (2001, 
PA = 0°, 45°, 90°), who inferred the presence of a supermas- 
sive black hole of 4 x 10 7 M@ by three-integral, axisymmet- 
ric modelling. Sil’chenko (1999) observed the circumnuclear 
region of NGC 1023 via two-dimensional spectroscopy, dis¬ 
entangling a kinematically and chemically decoupled stellar 
disc (with a radius of 80 pc). The stars of this nuclear disc 
have a mean age of 7 Gyr and are substantially younger than 
the rest of the galaxy, which is characterized by extremely 
red integrated colours (RC3) and without global star forma¬ 
tion (Pogge & Eskridge 1993). 


4 SURFACE PHOTOMETRY AND DATA 
REDUCTION 

4.1 Observations and reduction 

We observed NGC 1023 during December 27-28, 2000 at 
the 1-m Jacobus Kapteyn Telescope (JKT) located at the 
Roque de los Muchachos Observatory (ORM), La Palma. 
The detector consisted of a SITe2 device, with 2048 x 2048 
pixels and an image scale of 0 / . , 33 pixel -1 , giving an unvi¬ 
gnetted field of view of about 10' x 10'. The seeing during 
the observing run varied from 1" to l'/5. Tab. ^ gives the 
log of the photometric observations. The galaxy was imaged 
using the Harris B, V and I bandpasses. Deep observations 
were taken in all filters, reaching p/ ~ 22 mag arcsec -2 and 
/is ~ 24 mag arcsec -2 . 

The images were reduced using standard IRAF^J tasks. 
A run of 10 bias images was obtained each night, which were 
combined into a single bias frame and subtracted from the 
object frames. The images were flat-fielded using sky flats 
taken in all filters at the beginning and end of each observing 
night. The sky background level was removed by fitting a 

t IRAF is distributed by NOAO, which is operated by AURA Inc., 
under contract with the National Science Foundation 


Table 1. Log of the surface photometry observations 


Filter 

Date 

t exp 

00 

B 

27 Dec 2000 

3 x 900 


28 Dec 2000 

900 

V 

28 Dec 2000 

4 x 480 

I 

27 Dec 2000 

20 x 180 


28 Dec 2000 

10 x 180 


second order polynomial to the regions free of sources in the 
images. 

For the photometric calibration of the galaxies, stan¬ 
dard stars of known magnitudes were observed. The cali¬ 
bration constant includes corrections for atmospheric and 
Galactic extinction, and a colour term. No attempt was 
made to correct for internal extinction. The atmospheric 
extinction was taken from the differential aerosol extinction 
for the ORM (King 1985). The Galactic extinction in the 
B-band was taken from Schlegel et al. (1998). We used the 
Galactic extinction law from Cardelli et al. (1989) in order 
to get the extinction in the other filters. 

Fig. 0 shows the reduced image of the galaxy in the 
B-band and the B — I colour map image. 

4.2 Isophotal analysis and bar length 

We analysed the isophotal profiles of the galaxy by masking 
NGC 1023A and the foreground stars, then fitting ellipses 
using the IRAF task ELLIPSE. We hrst allowed the centers of 
the ellipses to vary, to test whether the optical disc is dis¬ 
turbed. Within the errors of the fits, we found no evidence of 
a varying center. Thus the inner disc of NGC 1023 has had 
enough time to settle if it had been disturbed in the tidal 
interaction suggested by the H I. The ellipse fits were there¬ 
fore repeated with the ellipse centers fixed; the resulting 
azimuthally averaged surface brightness, B — I colour, ellip- 
ticity and PA profiles are plotted in Fig. B. These are in good 
agreement with the surface photometry of Sanchez-Portal et 
al. (2000), although our surface photometry reaches fainter 
surface brightness. The PA and ellipticity profiles are similar 
for all bandpasses, suggesting that there is little, or uniform, 
obscuration, as required for the TW method. Furthermore, 
at large radii, the B—I profile (see Fig. |jj) is flat, which is fur¬ 
ther evidence of little obscuration. The only structure seen 
in the B — I profile is a small gradient in the region where 
the bulge and bar dominate. This gradient is in the sense of 
a bulge redder than the disc, which is most likely caused by 
an age difference between the bulge and disc stellar popula¬ 
tions. There is also a small colour gradient in the bar region, 
but the color variation is not patchy and is therefore unlikely 
to be caused by dust. We fitted the PA and inclination of the 
disc by averaging the /-band data at semi-major axis greater 
than 80”, finding PA = 80?2±0?5 and i = 66?4±1?2. Appli¬ 
cation of the TW method generally requires a PA measured 
to ~ 5° accuracy (Debattista & Williams 2002); we note 
that the value of the PA given in the RC3 is 87°. Our value 
of the PA is in very good agreement with that measured by 
Mollenhoff & Heidt (2001), who found PA = 79?15 in the 
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NGC 1023 B 



NGC 1023 



Figure 1. The B-band image (top) and B — I colour map (bot¬ 
tom) of NGC 1023. In the B-band image, the overplotted full line 
is the major axis (PA = 80?2), while the minor axis is indicated 
by the dot-dashed line. The offset slits, parallel to the major axis 
and shifted by 12" to the north and 16" to the south, are shown 
by dashed lines. The companion galaxy NGC 1023A can be seen 
to the east, just south of the major axis. The contour levels shows 
are /ib = 18 to 22 mag arcsec -2 , at unit intervals, with the outer¬ 
most contour being fig = 22 mag arcsec -2 . The B — I map shows 
that the strongest colour variations are associated with the bulge 
(which is red) and NGC 1023A (which is blue). 


J-band, 80?42 in the H- band, and 78?84 in the iF-band, for 
an average of 79? 5 ± 0?7. 

We measured the disc exponential scale-length, Rd, 
from our /-band surface photometry in the region outside 
the bar (R > 80”). We obtained Rd = 59”, corresponding 
to 2.9 kpc at our assumed distance. 

The most natural way to measure ag is with a Fourier 
analysis of the azimuthal luminosity profile (Ohta et al. 
1990; Aguerri et al. 2000), which we have applied here. We 
began by deprojecting the /-band image of the galaxy by 
a flux-conserving stretch along the minor axis by the factor 
1/cosi. NGC 1023 contains a significant spheroidal bulge: 
the bulge-disc model decompositions of Mollcnhoff & Heidt 
( 2001 ) for this galaxy have a bulge that accounts for about 
30% to 50% of the near-infrared light. Since we have not at¬ 
tempted to subtract the bulge, we expect that, in deprojec¬ 
tion, it appears as a structure elongated with the minor axis, 
dominating the center of the galaxy. We then decomposed 


the deprojected luminosity profile, /(/?,/>), where (/?,</) are 
plane polar coordinates in the galaxy frame, into a Fourier 
series: 

, . OO 

/(//,/>) = -3 3 — -(A m (_R) cos {m4>)+Bm{R) sin(m</))(3) 

m =1 

where the coefficients are defined by: 

1 [ 2 " 

Am(R) = — / I(R,(p) cos(m</)d/> (4) 

./n 


Bm(R) = — / /(//,/>) sin(m/>)d/> (5) 

n Jo 

The Fourier amplitude of the m-th component is defined as: 

( Aq(R)/ 2 m = 0 

I™(R) = { _ (6) 

[ \ZA 2 m (R) + B^R) ra/0 

Fig. shows the first (m = 1, 2, 3,4) relative Fourier ampli¬ 
tudes, Im/Io, of the /-band image as a function of R. The 
bar is evidenced by a strong m = 2 component. Fig. ^ also 
shows the bar/inter-bar intensity as a function of R. The bar 
intensity, h, is defined as the sum of the even Fourier com¬ 
ponents, Io + h, while the inter-bar intensity, Ub, is given 
by Io — I 2 (Ohta et al. 1990; Elmegreen & Elmegreen 1990; 
Aguerri et al. 2000). Ohta et al. (1990) arbitrarily defined 
ag as the outer radius for which h/Iib = 2 ; as Aguerri et al. 
( 2000 ) pointed out, a fixed value of h/Iib cannot account 
for the wide variety of bar luminosities present in galaxies. 
Instead, Aguerri et al. (2000) defined ag as the full width at 
half maximum (FWHM) of the curve of h/Iib■ Using this 
definition, but excluding the peak in h/Iib resulting from 
the deprojected bulge, we obtain ag = 69” ± 5”. The phase 
angle, 6 ^ 2 (/?) = tan _ 1 (A 2 (/?)//? 2 (//)), of the m = 2 compo¬ 
nent is shown in Fig. [| The bulge dominates in the inner 
20” (with 82 ~ 90°, as expected), with a transition from 
bulge to bar occurring in the range 20” ^ R 50”. In the 
region 50” ^ R ^ 70”, 82 has a minimum, implying the bar 
extends to somewhere here. However, 82 is not constant, be¬ 
cause of the influence of the bulge and the disc, and it cannot 
be used to determine ag. 

A independent measure of ag was obtained by depro¬ 
jecting the ellipses which best fit the isophotes, averaging 
over the three bands, and examining the deprojected ellip- 
ticity, e', as shown in Fig. [^. We find a sharp break in the 
slope of e at 69” ±3”, which we identify as the end of the bar. 
We repeated this experiment, deprojecting with different i 
and PA within the standard errors of these quantities. The 
errors in PA do not significantly change the semi-major axis 
at which the break occurs, while those in i lead to changes 
in the break semi-major axis of ~ 3”. Thus, our second es¬ 
timate of ag = 69” ± 4”. Fig. |f| plots the projected and 
deprojected ellipses, from which the change at semi-major 
axis ~ 70” can be discerned. 

Outside the inner 5”, the largest twist due to the bar 
of the isophotes away from the sky-plane major axis is only 
~ 8 °, which occurs at 40”. We measured the bar’s intrinsic 
angle from the major axis, i/bar, in the radial range 55” to 
70”, averaged over the 3 bands. Due to the contamination 
from the bulge, the deprojected isophotal ellipses are not 
constant in this range, so we used as the uncertainty in iph&T 
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Figure 2. Surface brightness (top-left), colour (top-right), ellipticity (center-left) and PA (center-right) radial profiles of the ellipse fits 
to the isophotes of NGC 1023. The full, dotted and dashed lines show results for /, V and B-bands respectively. The bottom two panels 
show the /—band ellipticity (left) and PA (right) for R > 80”, which we averaged to determine the disc’s inclination and PA. The solid 
lines represent the average values in this range: e = 0.595 (corresponding to i = 66?4) and PA = 80?2. 


the largest deviation from the average. Thus we obtained 
ipbn .1 = 102 ±2°, so that the bar is only some 12° from the 
minor axis, as can be seen in Fig. [I| This bar orientation, 
together with the large inclination and strong bulge, tend to 
obstruct our view of the bar. In this sense, NGC 1023 is not 
an ideal candidate for the TW method, and the expected 
values of V and X are likely to be small, as we show in 
0 § For this combination of inclination and i/>bar, Debattista 
(2002) suggests an accuracy of ~ 2? 5 in PA is required for 
accurate measurement of 


5 LONG-SLIT SPECTROSCOPY 


5.1 Observations and data reduction 

The spectroscopic observations of NGC 1023 were carried 
out at the ORM in La Palma, with the 3.6-m Telescopio 
Nazionale Galileo on February 01-02, 2001. The telescope 
mounted the Low Resolution Spectrograph (DOLORES). 
The HR-V grism No. 6 with 600 grooves mm -1 was used 
in combination with the 0”7 x 8ll slit and the thinned 
and back-illuminated Loral CCD with 2048 x 2048 pixels 
of 15 x 15 (im 2 . The wavelength range between 4697 and 
6840 A was covered with a reciprocal dispersion of 1.055 A 
pixel” 1 , which guarantees an adequate oversampling of the 
instrumental broadening function. Indeed, the instrumen¬ 
tal resolution, obtained by measuring the width of emission 
lines of a comparison spectrum after the wavelength calibra¬ 
tion, was 3.10 A (FWHM). This corresponds to an instru¬ 
mental dispersion a = 1.32 A (i.e., ~ 80 and ~ 60 km s” 1 
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Figure 3. Fourier analysis of the deprojected /-band light distribution of NGC 1023. The top-left panel shows the relative amplitudes of 
the m = 1 to m = 4 Fourier components, which reveals that m = 2 dominates. The top-right panel shows the phase angle of the m = 2 
Fourier component; inside R — 30", the system is clearly dominated by the spheroidal bulge, which deprojects into a feature along the 
minor axis. In the range 50” 7 R 5 70", 6 2 reaches a minimum of ~ 65°, but is not constant anywhere. The panel on the bottom-left 
shows the bar/inter-bar intensity ratio. The dashed line indicates one-half of the peak intensity (excluding the peak at R < 30” since 
this is an artifact of the bulge deprojection). The panel at bottom-right shows an alternative measure of ag, based on the deprojected 
ellipticities, e'\ the sharp break in e' at semi-major axis ~ 69” (indicated by the vertical dashed line) is identified as ag. In this panel, 
the upper, central and lower set of points correspond to a deprojection assuming i = 65?2, 66? 4 and 67?6 respectively. 


at the blue and red edges of the spectra, respectively). The 
angular sampling was 0”275 pixeR 1 . 

We obtained two 30-minutes spectra with the slit along 
the major axis (PA = 80°) and two offset spectra with the 
slit parallel to the major axis and shifted by 12” northward 
and 16” southward. The exposure time of each offset spec¬ 
trum was 60 minutes. During the observing run, we took 
spectra of the giant stars HR 2035 (G8III), HR 2429 (K1III), 
HR 2503 (K4III), HR 2701 (K0III), and HR 5370 (K3III) 
to use as templates in measuring the stellar kinematics. At 
the beginning of each exposure, the slit was positioned by 
acquiring with DOLORES a series of target images. A spec¬ 
trum of the comparison helium arc lamp was taken after 


each target exposure to allow an accurate wavelength cali¬ 
bration. The value of the seeing FWHM during the galaxy 
exposures ranged between O'!7 and I”0 as measured by fit¬ 
ting a two-dimensional Gaussian to the guide star. 

All the spectra were bias subtracted, Hatfield corrected, 
cleaned of cosmic rays, and wavelength calibrated using 
standard MIDAS0 routines. The bias level was determined 
from the bias frames obtained during the observing nights 
to check the CCD status. The Hatfield correction was per¬ 


il MIDAS is developed and maintained by the European Southern 
Observatory 
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12" North Major axis 16" South 



Figure 5. Stellar kinematics of NGC 1023. Each column shows the radial profiles of the four Gauss-Hermite moments (V, cr, /13 and / 1 . 4 ) 
we measured along a given slit position. In the top panels the horizontal dashed line indicates the heliocentric systemic velocity of the 
NGC 1023. All velocities are plotted as observed, without correcting for the inclination. The full data are given in Tab. I 2 I 


formed by means of both quartz lamp and twilight sky spec¬ 
tra, which were normalized and divided into all the spectra, 
to correct for pixel-to-pixel sensitivity variations and large- 
scale illumination patterns due to slit vignetting. Cosmic 
rays were identified by comparing the counts in each pixel 
with the local mean and standard deviation (as obtained 
from Poisson statistics by taking into account the gain and 
read-out noise of the detector) and then corrected by in¬ 
terpolating over. The residual cosmic rays were corrected 
by manually editing the spectra. The wavelength calibra¬ 
tion was performed by means of the MIDAS package XLONG. 
Each spectrum was rebinned using the wavelength solu¬ 
tion obtained from the corresponding arc-lamp spectrum. 
We checked that the wavelength rebinning had been done 
properly by measuring the difference between the measured 
and predicted wavelengths (Osterbrock et al. 1996) for the 
brightest night-sky emission lines in the observed spectral 
range. The resulting accuracy in the wavelength calibration 
is ~ 1 km s -1 . The major-axis spectra were co-added using 


the center of the stellar continuum as reference. In the re¬ 
sulting spectra, the contribution of the sky was determined 
by interpolating along the outermost 20” at the two edges 
of the slit, where the galaxy light was negligible, and then 
subtracted. 

5.2 Stellar kinematics 

We measured the stellar kinematics from the galaxy ab¬ 
sorption features present in the wavelength range running 
from 4697 A to 5570 A and centered on the Mg line triplet 
(AA 5164, 5173, 5184 A). The spectra were rebinned along the 
spatial direction to obtain a nearly constant signal-to-noise 
ratio larger than 20 per resolution element (with a peak of 
50 in the innermost regions of the major axis spectrum). 
At the outermost radii the signal-to-noise ratio decreases to 
~ 10. The galaxy continuum was removed row-by-row by 
fitting a fourth to sixth order polynomial. 

We used the Fourier Correlation Quotient method 
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Figure 4. Best-fitting ellipses to the /-band isophotes (a) before 
and (b) after deprojection to the NGC 1023 disc plane (i = 06?4 
and PA= 80? 2). The system has been rotated so that the disc 
major axis is horizontal. Bold-line ellipses have deprojected semi¬ 
major axes of 00". 70" and 80"; the straight lines indicate the 
major axes of these 3 ellipses. 


(FCQ, Bender 1990) following the prescriptions of Bender, 
Saglia & Gerhard (1994) and adopting HR 2035 as kinemat- 
ical template. This allowed us to derive, for each spectrum, 
the line-of-sight velocity distribution (LOSVD) along the slit 
and to measure its moments, namely the radial velocity v, 
the velocity dispersion a and the values of the coefficients 
hi and h 4 . At each radius, they have been derived by fitting 
the LOSVD with a Gaussian plus third- and fourth-order 
Gauss-Hermite polynomials H 3 and H 4 , which describe the 
asymmetric and symmetric deviations of the LOSVD from 
a pure Gaussian profile (van der Marel & Franx 1993; Ger¬ 
hard 1993). We derived errors on the LOSVD moments from 
photon statistics and CCD read-out noise, calibrating them 
by Monte Carlo simulations as done by Bender et al. (1994) 
and Bower et al. (2001), who have recently derived the stel¬ 
lar kinematics of this galaxy along three axes (different from 
ours). In general, our errors are in the range of 10-20 km s _1 
for v and a, and of 0.03-0.1 for hi and hi, becoming larger 
for signal-to-noise ratios lower than 20. These errors do not 
take into account possible systematic effects due to any tem- 
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In 00 


Figure 6. Galaxy spectra (thick lines), obtained by collapsing 
along the spatial direction, of the major axis and the two offset 
positions. They are compared with the template spectra (thin 
line) convolved with their corresponding LOSVDs. Galaxy and 
template spectra have been continuum subtracted and tapered at 
the ends with a cosine bell function. 

plate mismatch. The measured stellar kinematics along the 
major and the offset axes are reported in Tab. ^ and plotted 
in Fig. §. 


6 PATTERN SPEED 

To compute the mean position of stars, X, along the slits, we 
extracted profiles from the B, V and /-band surface photom¬ 
etry along the positions of the slits. The profiles in V-band 
(within which the wavelength range used in the spectroscopy 
falls) match very well the profiles obtained by collapsing the 
spectra along the wavelength direction, confirming that the 
slits were placed as intended. We used the broad band pro¬ 
files to compute X because these are less noisy than the 
spectral profiles, particularly at large radii. We therefore 
computed three values of X at each slit position, Xb , Xy 
and Xi. Each value of X was computed by Monte Carlo 
simulation, with photon, read-out and sky noise to com¬ 
pute the errors. Another advantage to using the broad band 
surface photometry is that it allows us to test for contami¬ 
nation by NGC 1023A to the integrals. Formally, the inte¬ 
grals in Eqn. [I] are over —oo < X < oo, but can be lim¬ 
ited to —X m ax < X < Xmax if Xmax has reached the ax- 
isymmetric part of the disc; still larger values of X max add 
noise only. In the case of NGC 1023, the presence of NGC 
1023A introduces some additional light at negative X. The 
most severely contaminated slit is the one at 16"S, where X 
continues changing as ever larger values of X max are used. 
We therefore found it useful to compare the 16"S profiles 
to the profiles at 16"N. For these latter profiles, there is 
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Table 2. Stellar kinematics of NGC 1023. This data is plotted in Fig. 
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25.89 

497.4 

± 

16.8 

164.7 ± 18.2 

0.068 

± 

0.091 

- 0.044 

± 

0.106 




Major Axis 





31.20 

475.4 

± 

18.1 

137.7 ± 

15.5 

0.063 

± 

0.136 

- 0.072 

± 

0.128 











38.25 

438.6 

± 

15.5 

115.9 ± 11.8 

0.082 

± 

0.127 

- 0.067 

± 

0.107 

- 108.76 

852.3 ± 

42.5 

95.8 ± 

27.9 

- 0.090 ± 

0.232 

- 0.365 

± 

0.200 

48.38 

430.4 

± 

11.4 

82 .- 

4 ± 9.9 

- 0.020 

± 

0.105 

- 0.073 

± 

0.071 

- 82.56 

834.8 ± 

28.5 

98.1 ± 

19.2 

- 0.109 ± 

0.200 

- 0.171 

± 

0.271 

64.33 

403.5 

± 

19.9 

85.2 ± 13.6 

0.119 

± 

0.165 

- 0.148 

± 

0.130 

- 59.33 

791.9 ± 

24.5 

114.1 ± 

20.2 

- 0.040 ± 

0.194 

- 0.060 

± 

0.103 

89.29 

407.5 

± 

28.0 

92.6 ± 19.0 

0.148 

± 

0.199 

- 0.161 

± 

0.327 

- 45.82 

808.1 ± 

18.5 

107.3 ± 

15.4 

- 0.091 ± 

0.150 

- 0.085 

± 

0.118 

112.90 

392.2 

± 

45.5 

84.7 ± 30.8 

0.311 

± 

0.277 

- 0.172 

± 

0.356 

- 36.62 

783.6 ± 

16.4 

123.6 ± 

13.2 

- 0.040 ± 

0.151 

- 0.053 

± 

0.123 














- 29.99 

738.7 ± 

19.4 

154.1 ± 

16.2 

- 0.063 ± 

0.112 

- 0.072 

± 

0.127 
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North 

Offset 






- 24.93 

725.0 ± 

15.7 

168.3 ± 

15.2 

0.046 ± 

0.078 

- 0.082 

± 

0.099 














- 20.98 

720.3 ± 

17.7 

164.9 ± 

18.3 

- 0.019 ± 

0.066 

- 0.003 

± 

0.104 

- 78.35 

814.6 

± 

41.8 

121.2 

± 

50.2 

- 0.057 

± 

0.233 

- 0.073 

± 

0.231 

- 17.83 

717.0 ± 

17.0 

164.8 ± 

14.6 

- 0.023 ± 

0.086 

- 0.104 

± 

0.114 

- 56.44 

809.1 

± 

32.8 

109.5 

± 

28.5 

0.024 

± 

0.211 

- 0.107 

± 

0.170 

- 15.37 

716.4 ± 

13.7 

166.8 ± 

14.1 

0.047 ± 

0.064 

- 0.057 

± 

0.087 

- 43.07 

760.6 

± 

30.4 

111.4 

± 

22.5 

- 0.109 

± 

0.217 

- 0.107 

± 

0.145 

- 13.32 

717.4 ± 

14.4 

168.9 ± 

16.1 

0.028 ± 

0.059 

- 0.012 

± 

0.086 

- 34.80 

743.5 

± 

28.5 

126.2 

± 

18.8 

0.046 

± 

0.140 

- 0.097 

± 

0.194 

- 11.53 

713.1 ± 

13.1 

160.8 ± 

12.9 

- 0.041 ± 

0.057 

- 0.025 

± 

0.083 

- 28.68 

714.1 

± 

27.6 

125.0 

± 

17.1 

0.020 

± 

0.137 

- 0.162 

± 

0.194 

- 10.02 

703.1 ± 

12.4 

191.4 ± 

12.3 

0.014 ± 

0.043 

- 0.029 

± 

0.048 

- 23.61 

709.2 

± 

27.1 

118.7 

± 

21.7 

- 0.014 

± 

0.168 

- 0.079 

± 

0.206 

- 8.79 

685.3 ± 

12.2 

192.1 ± 

10.5 

0.007 ± 

0.043 

- 0.071 

± 

0.051 

- 19.23 

669.9 

± 

23.6 

154.7 

± 

17.6 

0.078 

± 

0.122 

- 0.126 

± 

0.145 

- 7.69 

693.1 ± 

13.9 

189.1 ± 

14.0 

0.002 ± 

0.048 

- 0.030 

± 

0.056 

- 15.39 

682.0 

± 

17.3 

146.6 

± 

14.3 

0.019 

± 

0.088 

- 0.100 

± 

0.129 

- 6.74 

695.1 ± 

12.2 

191.3 ± 

11.7 

- 0.010 ± 

0.041 

- 0.024 

± 

0.045 

- 12.11 

653.2 

± 

15.9 

143.7 

± 

13.6 

0.010 

± 

0.081 

- 0.078 

± 

0.119 

- 5.91 

691.8 ± 

11.0 

173.2 ± 

12.2 

- 0.019 ± 

0.043 

- 0.001 

± 

0.058 

- 9.22 

650.4 

± 

19.0 

149.4 

± 

16.8 

- 0.015 

± 

0.086 

- 0.047 

± 

0.118 

- 5.08 

668.9 ± 

12.9 

191.4 ± 

12.8 

- 0.024 ± 

0.043 

- 0.017 

± 

0.048 

- 6.49 

636.1 

± 

17.5 

171.1 

± 

17.5 

0.081 

± 

0.078 

- 0.032 

± 

0.084 

- 4.40 

670.7 ± 

13.0 

181.4 ± 

14.9 

0.013 ± 

0.047 

0.004 

± 

0.059 

- 3.88 

626.6 

± 

17.8 

186.2 

± 

19.7 

0.003 

± 

0.066 

0.000 

± 

0.074 

- 3.85 

660.2 ± 

12.6 

200.7 ± 

12.1 

- 0.005 ± 

0.042 

- 0.022 

± 

0.042 

- 1.40 

602.2 

± 

14.5 

176.6 

± 

14.3 

0.033 

± 

0.057 

- 0.028 

± 

0.069 

- 3.30 

665.0 ± 

10.6 

200.2 ± 9.9 

- 0.032 ± 

0.035 

- 0.013 

± 

0.035 

1.07 

592.4 

± 

14.7 

159.9 

± 

13.8 

0.036 

± 

0.064 

- 0.034 

± 

0.081 

- 2.75 

669.6 ± 

11.5 

196.0 ± 

10.3 

- 0.022 ± 

0.039 

- 0.036 

± 

0.040 

3.53 

585.8 

± 

16.9 

167.8 

± 

17.3 

0.029 

± 

0.069 

- 0.018 

± 

0.083 

- 2.34 

681.8 ± 

12.7 

194.3 ± 

12.9 

- 0.028 ± 

0.041 

- 0.003 

± 

0.045 

6.13 

539.7 

± 

16.3 

177.7 

± 

15.6 

0.069 

± 

0.067 

- 0.056 

± 

0.080 

- 2.07 

674.1 ± 

11.3 

198.6 ± 

12.5 

0.004 ± 

0.037 

0.019 

± 

0.041 

9.02 

522.3 

± 

19.0 

152.1 

± 

17.1 

- 0.063 

± 

0.085 

- 0.043 

± 

0.112 

- 1.79 

672.3 ± 

11.8 

199.6 ± 

11.1 

- 0.040 ± 

0.039 

- 0.019 

± 

0.039 

12.29 

542.2 

± 

15.9 

162.3 

± 

16.6 

0.036 

± 

0.067 

0.004 

± 

0.082 

- 1.52 

674.0 ± 

10.3 

203.1 ± 

10.7 

- 0.030 ± 

0.033 

0.017 

± 

0.035 

16.12 

547.4 

± 

18.2 

176.4 

± 

18.7 

0.036 

± 

0.071 

- 0.030 

± 

0.084 

- 1.24 

683.3 ± 

11.3 

208.2 ± 

10.2 

- 0.027 ± 

0.038 

- 0.029 

± 

0.035 

20.76 

507.8 

± 

21.4 

131.6 

± 

14.1 

0.046 

± 

0.141 

- 0.122 

± 

0.154 

- 0.97 

681.2 ± 

11.2 

202.6 ± 

10.5 

- 0.030 ± 

0.037 

- 0.016 

± 

0.036 

26.72 

483.6 

± 

29.1 

138.1 

± 

20.0 

0.048 

± 

0.136 

- 0.122 

± 

0.162 

- 0.69 

661.3 ± 

12.2 

224.3 ± 

11.1 

- 0.008 ± 

0.042 

- 0.029 

± 

0.037 

34.99 

446.8 

± 

24.0 

99.0 

± 

14.0 

0.058 

± 

0.213 

- 0.191 

± 

0.132 

- 0.42 

648.1 ± 

10.5 

216.1 ± 

10.5 

- 0.024 ± 

0.034 

0.005 

± 

0.033 

47.61 

416.4 

± 

57.6 

117.4 

± 

46.8 

0.158 

± 

0.236 

- 0.122 

± 

0.192 

- 0.14 

615.9 ± 

23.2 

247.0 ± 

32.0 

0.020 ± 

0.131 

0.041 

± 

0.113 

66.86 

414.1 

± 

27.1 

68.4 

± 

20.1 

- 0.058 

± 

0.255 

- 0.161 

± 

0.119 

0.13 

609.5 ± 

25.3 

236.3 ± 

29.6 

0.064 ± 

0.100 

- 0.002 

± 

0.076 

96.37 

424.1 

± 

32.3 

55.7 

± 

31.0 

0.001 

± 

0.230 

- 0.167 

± 

0.165 

0.41 

558.6 ± 9.2 

215.1 ± 

10.2 

0.124 ± 

0.031 

- 0.002 

± 

0.036 














0.68 

550.7 ± 9.4 

215.5 ± 9.7 

0.082 ± 

0.032 

- 0.001 

± 

0.033 
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South 

Offset 






0.96 

534.0 ± 

10.6 

211.9 ± 

10.7 

0.062 ± 

0.037 

- 0.011 

± 

0.037 














1.23 

532.8 ± 

10.1 

200.7 ± 

10.8 

0.066 ± 

0.036 

- 0.011 

± 

0.040 

- 96.43 

834.8 

± 

38.4 

79.1 

± 

25.4 

- 0.075 

± 

0.250 

- 0.270 

± 

0.179 

1.51 

536.2 ± 

10.1 

205.6 ± 

11.3 

0.083 ± 

0.035 

0.003 

± 

0.039 

- 75.58 

800.9 

± 

37.1 

77.9 

± 

21.6 

0.006 

± 

0.199 

- 0.311 

± 

0.223 

1.78 

527.2 ± 

11.7 

211.4 ± 

10.9 

0.019 ± 

0.040 

- 0.028 

± 

0.038 

- 53.66 

793.2 

± 

29.8 

94.7 

± 

26.1 

- 0.029 

± 

0.206 

- 0.079 

± 

0.140 

2.06 

527.3 ± 

15.7 

211.5 ± 

16.4 

0.020 ± 

0.053 

- 0.006 

± 

0.052 

- 38.07 

779.7 

± 

28.9 

88.7 

± 

18.8 

0.029 

± 

0.224 

- 0.109 

± 

0.128 

2.33 

529.1 ± 

12.2 

213.2 ± 

12.0 

0.052 ± 

0.042 

- 0.012 

± 

0.040 

- 28.42 

786.6 

± 

28.3 

100.5 

± 

16.7 

0.020 

± 

0.218 

- 0.107 

± 

0.112 

2.61 

534.5 ± 

13.6 

206.3 ± 

16.4 

0.116 ± 

0.048 

0.019 

± 

0.056 

- 21.04 

711.6 

± 

30.2 

152.8 

± 

27.2 

- 0.053 

± 

0.178 

- 0.086 

± 

0.135 

3.01 

552.2 ± 

12.0 

207.6 ± 

15.6 

0.097 ± 

0.039 

0.066 

± 

0.048 

- 15.18 

704.1 

± 

28.8 

134.2 

± 

27.2 

- 0.052 

± 

0.186 

0.009 

± 

0.128 

3.56 

550.0 ± 

11.7 

210.8 ± 

12.0 

0.040 ± 

0.040 

- 0.004 

± 

0.040 

- 10.40 

647.1 

± 

27.0 

154.9 

± 

21.1 

- 0.094 

± 

0.184 

- 0.152 

± 

0.145 

4.11 

549.5 ± 

11.3 

211.7 ± 

11.4 

0.058 ± 

0.039 

- 0.009 

± 

0.039 

- 6.15 

654.5 

± 

21.9 

145.5 

± 

20.2 

- 0.029 

± 

0.140 

- 0.040 

± 

0.132 

4.66 

533.6 ± 

11.7 

194.4 ± 

13.4 

0.046 ± 

0.041 

0.009 

± 

0.047 

- 2.18 

633.8 

± 

30.8 

164.8 

± 

24.3 

- 0.104 

± 

0.158 

- 0.134 

± 

0.162 

5.35 

538.3 ± 9.0 

196.4 ± 

10.7 

0.111 ± 

0.036 

- 0.038 

± 

0.045 

1.52 

606.4 

± 

23.7 

163.1 

± 

18.2 

0.150 

± 

0.131 

- 0.186 

± 

0.154 

6.17 

531.4 ± 

13.4 

204.0 ± 

14.8 

- 0.017 ± 

0.043 

0.017 

± 

0.046 

5.23 

594.6 

± 

18.0 

140.3 

± 

16.0 

- 0.056 

± 

0.140 

- 0.057 

± 

0.113 

7.00 

534.5 ± 

12.6 

191.2 ± 

12.8 

0.013 ± 

0.043 

- 0.022 

± 

0.049 

9.06 

594.3 

± 

24.9 

162.5 

± 

24.2 

- 0.001 

± 

0.101 

- 0.039 

± 

0.144 

7.95 

537.8 ± 

14.2 

189.2 ± 

14.5 

- 0.024 ± 

0.048 

- 0.021 

± 

0.056 

13.05 

591.2 

± 

21.2 

150.5 

± 

19.5 

0.020 

± 

0.133 

- 0.064 

± 

0.142 

9.06 

525.0 ± 

14.3 

203.0 ± 

13.4 

0.006 ± 

0.048 

- 0.035 

± 

0.047 

17.43 

548.2 

± 

25.6 

133.7 

± 

21.7 

0.017 

± 

0.189 

- 0.064 

± 

0.127 

10.29 

518.0 ± 

13.7 

188.2 ± 

14.9 

0.007 ± 

0.047 

- 0.012 

± 

0.056 

22.37 

516.1 

± 

23.3 

128.5 

± 

19.2 

0.072 

± 

0.164 

- 0.064 

± 

0.128 

11.80 

533.2 ± 

13.9 

186.6 ± 

15.0 

0.043 ± 

0.051 

- 0.034 

± 

0.063 

27.83 

511.8 

± 

24.2 

129.7 

± 

19.3 

0.061 

± 

0.177 

- 0.096 

± 

0.134 

13.58 

559.2 ± 

12.8 

186.9 ± 

13.4 

0.103 ± 

0.055 

- 0.087 

± 

0.072 

34.42 

489.1 

± 

37.7 

141.7 

± 

34.6 

0.296 

± 

0.191 

0.040 

± 

0.204 

15.77 

531.0 ± 

14.1 

173.2 ± 

16.3 

- 0.004 ± 

0.053 

0.014 

± 

0.075 

43.89 

458.7 

± 

35.7 

108.3 

± 

27.8 

0.135 

± 

0.250 

- 0.075 

± 

0.147 

18.50 

509.1 ± 

15.1 

160.1 ± 

14.4 

0.013 ± 

0.076 

- 0.071 

± 

0.104 

59.92 

438.6 

± 

36.8 

85.1 

± 

25.3 

- 0.001 

± 

0.198 

- 0.184 

± 

0.179 

21.79 

502.6 ± 

17.8 

151.0 ± 

16.8 

0.123 ± 

0.113 

- 0.053 

± 

0.106 

86.58 

489.2 

± 

50.1 

89.8 

± 

44.6 

0.183 

± 

0.212 

- 0.124 

± 

0.235 


no further signal to X beyond X max = 60". The compari¬ 
son also showed that there is little or no contamination to 
the 16”S slit inside X max ~ 55”. Furthermore, since NGC 
1023A is quite blue, the contamination decreases in going 
from B to V, and especially to / band. Therefore, to mea¬ 
sure X at 16”S, we averaged the Monte Carlo results be¬ 
tween 62” < Xmax < 75” (75” < R < 85”, where the disc is 
roughly axisymmetric) and then averaged over each band. 
In the remaining slits, contamination by NGC 1023A in the 
major axis begins at R > 100” in all bands and is easily 
avoided, while the 12”N slit has no obvious contamination. 
We therefore also averaged them over 75” < R < 85”. Be¬ 
cause of the color gradient in the inner ~ 70” (Fig. |2|), the 
value of X in the offset slits changes with the color band 
used, with \Xi\ > \Xy\ \Xb\- Therefore, by averaging over 
bands and then taking the largest discrepancy from the av¬ 
erage as our error value, we are confident that we take into 
account any difference between the narrow wavelength range 
of the spectroscopy and the broad one of the surface pho¬ 
tometry. The values thus obtained are given in Tab. [i| 

If we have successfully avoided contamination by NGC 
1023A, then we expect that, since NGC 1023 is nearly point- 


symmetric about its center, the profile extracted from the 
surface photometry at Y should have a value of X roughly 
the negative of value for the profile extracted at —Y. For a 
profile extracted from the surface photometry at 16”N, we 
found X = — 2”31 ± 0”15 averaged over the 3 bands, which 
is within 2cr of our value for —A at 16”S. 

To measure the luminosity-weighted line-of-sight stel¬ 
lar velocity, V, for each slit position, we collapsed each bi- 
dimensional spectrum along its spatial direction to obtain a 
one-dimensional spectrum. The resulting spectra have been 
analysed with the FCQ method using HR 2035 as template 
star. V is the radial velocity derived from the LOSVD of 
the one-dimensional spectra (see Fig. |ij). For each slit po¬ 
sition, the uncertainties on V were estimated by measur¬ 
ing it in different wavelength ranges between 4910 A and 
5560 A (but always including the Mg triplet). The V val¬ 
ues we derived along each slit are given in Tab. 0 We in¬ 
tegrated the spectral rows at .Y m „ = 140”, 100" and 80” 
in the major-axis, northern and southern offset spectra, re¬ 
spectively. For the major-axis and northern offset spectra we 
were limited by the noise, which becomes dominant at larger 
radii ( S/N < 3), while for the southern offset spectrum we 
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Figure 7. The pattern speed measurement for NGC 1023. The bottom 3 panels show the T-band profiles for the three slits (12”N, 
major axis and 16”S from left to right) with X indicated by the dashed vertical lines. The horizontal dotted lines indicate [iy = 19 mag 
arcsec -2 and fiy = 21 mag arcsec -2 . The 3 panels on the right indicate the LOSVD for the three slits (12"N, major axis and 1 6”S from 
top to bottom), with V indicated by the dashed lines. The LOSVDs are constructed from the Gauss-Hermite moments; the 3 LOSVDs in 
each of these panels were obtained from different wavelength ranges (here, each is renormalized to the same maximum value, for clarity). 
The central panel plots V against X (solid circles); the slope of the best fit straight line (indicated by the solid line) is Up sin i = 4.7±1.7 
km s —1 arcsec -1 . The open triangles and dotted line indicates the results for using just the /-band to measure X, which gives the largest 
values of \X\ and therefore the smallest slope, sin i = 4.0 ± 1.4 km s —1 arcsec -1 . 


decided to also try X max = 60”, to check for contamination 
by NGC 1023A. We found that different choices of X max do 
not affect our results for V. The values of V we derived by 
integrating both offset spectra out to = 140” (after 

removing by linear interpolation the contribution of fore¬ 
ground stars) are close to those in Tab. although their 
errors increase by a factor 1.6 and 2 along the north and 
south slit, respectively. Therefore, within the errors obtained 
from the varying wavelength range, there is no evidence of 
contamination by NGC 1023A in V for the southern offset 
slit. 

By performing our analysis with all the remaining tem- 


Table 3. The values of V and X for the 3 slits. 


Offset 

X 

V 


(”) 

(km s -1 ) 

12" N 

-1.82 ±0.13 

619.1 + 5.4 

Major axis 

+0.05 ±0.06 

614.5 + 2.2 

16” S 

+1.90 + 0.25 

603.7 + 3.0 


plate spectra, we found that V values are not affected by 
template mismatch (at least for our set of template stars). 
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Table 4. Heliocentric systemic velocity based on optical mea¬ 
surements. 


Vsys 

Source 

(km s -1 ) 


734 ± 41 

Mayall &; de Vaucouleurs (1962) 

617 ±19 

Tonry & Davies (1981) 

597 ± 9 

Schechter (1983) 

608 ± 9 

Schechter (1983) 

615 ± 20 

Dressier & Sandage (1983) 

620 ± 50 

Capaccioli et al. (1986) 

592 ± 15 

Simien & Prugniel (1997) 

606 ±5 

Weighted average of literature values 

613 ±3 

This work 


This is in agreement with the results of Bender et al. (1994) 
on the determination of the radial velocities with FCQ. 

Using the data of Tab. in Fig. (?] we plot V versus X, 
fitting a straight line (using subroutine fitexy in Numer¬ 
ical Recipes). The slope of this line is Uj,sini = 4.7 ± 1.7 
km s -1 arcsec -1 , which gives Q. p = 5.1±1.8 km s -1 arcsec -1 
(corresponding to D p = 104±34 km s -1 kpc -1 at 10.2 Mpc). 
We have also fitted a “slowest-bar” line, using the values of 
Xr (which are largest, and therefore give a shallower slope) 
in the 12"N and major-axis, and adopting the value of Xj 
from 16"N, which gives a slightly larger \X\ than the 16"S 
profile. The slope of this line, now ignoring the uncertainties 
in X, is Up sin i = 4.0 ±1.4 km s -1 arcsec -1 (U p = 4.3 ±1.5 
km s' 1 arcsec -1 ). 

To then determine 1Z, we need to measure Dl, which 
we approximate by the co-rotation radius, R c , from the ax- 
isymmetric approximation. Debattista (1998) found that, 
for strong bars in IV-body simulations, the approximation 
Dl — R c involves an error of ~ 5%. We therefore measured 
the rotation curve, V C {R), of NGC 1023. This is a two step 
process: first we fitted tilted rings (with constant PA and 
inclination fixed to the value obtained from the surface pho¬ 
tometry) to the data from all our slits, interpolating with 
cubic splines where necessary. This allowed us to measure 
the heliocentric systemic velocity, V sys = 613 ± 3 km s -1 , 
which is compared to the values based on optical data avail¬ 
able in the literature in Tab. [|. By folding the major axis 
slit (which reaches furthest out in R) about the origin, af¬ 
ter subtracting out V sys , we obtained the stellar streaming 
velocities, V*(R), shown in Fig. ||. 

We tested V*(R) by rotating it into the orientations of 
other slit observations available in the literature. The slits 
we compared with are those at PA = 87° (Simien & Prugniel 
1997; Neistein et al. 1999), PA = 90° (Bower et al. 2001) and 
PA = 177° (Simien & Prugniel 1997). Excluding the inner 
2 " , to avoid the effects of different seeing and/or slit widths, 
we find \ 2 ^ 1 in all cases except for the PA = 87° slit of 
Neistein et al. (1999), for which we obtained x 2 = 2.5 (for 
27 points). Since the same PA is also covered by a slit from 
Simien & Prugniel (1997), which we match well ( x 2 = 0-8 for 
36 points), we conclude that the stellar streaming velocity 
profile we obtained is reasonable and in good agreement with 
previous data. 

We then corrected for the asymmetric drift using the 
velocity dispersion data, to recover the true circular velocity, 



Figure 8 . The stellar velocity streaming curve, V*, derived by 
subtracting Vsys , folding and deprojecting the observed velocities 
on the major axis. The data from the major axis are shown in dia¬ 
monds (approaching side) and stars (receding side). The solid line 
and the dashed lines flanking it are our fitted spline together with 
the la error interval. The asymmetric drift terms (V/ 2 — V* 2 ) 
are shown by the dot-dashed, dotted, and dash triple-dotted lines 
(a = 0.0, 0.5, and 1.0 respectively) in the lower-right corner at 
R > 70". The corresponding circular rotation curves, V c , are 
shown in the bold lines. 


V C (R). Following Eqn. 4-33 in Binney & Tremaine (1987), 
and assuming 8 { vrv z ) / dz ~ 0, we can write the asymmetric 
drift equation as: 


T t/2 2 

V c — V* = -(J R 


din P , dIn a 2 R f 
9 In R a In R \ 



(7) 


The observed velocity dispersion, at a point ( R , </>) in the 
galaxy frame, is given by: 

of os = isin 2 ?. [(cr|± al + 2al cot 2 i) - (a% - al) cos 2<j>\ 
= a(R, <j>) + b(R, (/>) cos 2cj>. (8) 


where <j> = 0 is taken to be the disc’s major axis. When the 
disc is axisymmetric, a and b can only depend on R, in which 
case a and b are simply the first 2 even moments in a Fourier 
expansion of the line-of-sight velocity dispersion. For a disc 
with a flat rotation curve, ct^/ctr = 1/2, which then allows 
us to obtain from b. However, the profile of V * we found 
for NGC 1023 is not obviously flat. (Furthermore, b is more 
sensitive to noisy data than is a.) In the general case, the 
asymmetric drift equation can be re-written in terms of both 
a and b, as follows: 


R 


sin 2 i(l + a 2 cot 2 i) Y Rd 
2 b 


a — b d 

-dR [a ~ b) 


(9) 


where we have assumed that the disc density is exponen¬ 
tial with scale-length Rd and a = tj z /(jR. Using the velocity 
dispersion data from our 3 slits and the near minor-axis 
data of Simien & Prugniel (1997), we measured a — b and 
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R c) 


Figure 9. The average value, a , (top panel) and the cos 2 <f> com¬ 
ponent, h. (bottom panel) of the squared observed velocity dis¬ 
persions. The top panel also indicates a — b. We fit straight lines 
to a — b and b at R > 70”, as indicated. 

b. We fitted splines to interpolate between data points, and 
estimate errors from Monte Carlo simulation with Gaussian 
noise added from the error estimate on each point. The re¬ 
sulting profiles of a and b are shown in Fig. Our procedure 
for measuring V*(R), a and b is valid only in the axisym- 
metric region outside the bar. Therefore we only apply the 
asymmetric drift correction to R > 70”; for this region, we 
fit straight lines to a — b = —141 (±61) R + 41884 (±562) 
km 2 s -2 and to b = -319 (±39)1?. ± 8706 (±350) km 2 s -2 . 
With these values, we then correct the rotation curve, as¬ 
suming 0.0 < a < 1.0. (This range of a is a reasonable one: 
in the Solar neighborhood, Dehnen & Binney (1998) found 
a ~ 0.541q'q 2 , while in NGC 488, Gerssen et al. (1997) mea¬ 
sured a = 0.70 ± 0.19.) The resulting rotation curve is more 
or less flat, with 259 ^ 278 km s Adding in 

quadrature a velocity uncertainty of 30 km s _1 , we obtain 
V C! flat = 270 ± 31 km s _1 . (For comparison, Neistein et al. 
1999 found V c = 250 ± 17 km s -1 ; we have checked that, 
using a method similar to their Eqn. 8, gives similar values 
for V c jiat to those reported above.) If we extrapolate this 
value to smaller radii (which, at worst, over-estimates R c 
and, correspondingly, 1 Z , unless the rotation curve is drop¬ 
ping between 40” ^ R < 70”, which seems unlikely from 
Fig. ||) we find R c = V c ,flat/Q P = 53 Th! 2 ®”, where the 
first error is due to the uncertainties in V c jut and the sec¬ 
ond error reflects uncertainties in Q. p . We obtain, therefore, 
77 = 0.77±0.10+q2o (where the first error now also includes 
the uncertainty in an). For our “slowest-bar” estimate of the 
slope, we find 77 = 0.91 ± 0.121 q 23 . We have found, there¬ 
fore, that the bar in NGC 1023 is consistent with being a 
fast bar. 

7 DISCUSSION AND CONCLUSIONS 

We have found that the bar in NGC 1023 is fast, as are all 
bars which have been measured to date. Debattista & Sell- 


wood (1998, 2000) showed that fast bars can persist only if 
the disc is maximal. Following Ostriker & Peebles (1973), 
it is sometimes thought that the unbarred (SA) galaxies are 
stabilized by massive DM halos. However, massive DM halos 
are not necessary for stabilizing discs; a rapidly rising rota¬ 
tion curve in the inner disc, such as when a massive bulge is 
present, is also able to inhibit bar formation (Toomre 1981; 
Scllwood & Evans 2001). Debattista & Sellwood (1998) ar¬ 
gued that unbarred HSB galaxies must also be maximal for, 
if HSB disc galaxies form a continuum of DM halo masses 
spanning massive DM halo-stabilised SA galaxies to maxi¬ 
mal SB galaxies, then slow bars must also be found in the 
intermediate range of halo masses. If we seek to avoid inter¬ 
mediate halo masses and slow bars by postulating (for what¬ 
ever reason) a bimodal DM halo mass distribution for HSB 
galaxies, then we are left with the possibility that tidal in¬ 
teractions can still form bars, which would be slow (Noguchi 
1987; Salo 1991; Miwa & Noguchi 1998). Thus Debattista & 
Scllwood (2000) concluded that the absence of slow bars re¬ 
quires that all HSB disc galaxies are maximal. However, it 
is possible that no such slow bars have been found because 
of an observational bias against SB systems with evidence 
of tidal interactions. 

We have chosen to study NGC 1023, in part, because it 
shows signs of a weak interaction in its past, without being 
at present significantly perturbed. The fast bar we found 
indicates that NGC 1023 has a maximal disc. If SA galaxies 
are stabilized by massive halos, we should find slow bars in 
that fraction of SB galaxies in which the bar formed through 
the interaction. While it is not possible to reach a general 
conclusion on the DM content of SA galaxies based on our 
measurement for a single galaxy, a large enough sample of 
similar SB galaxies with mild interactions in the past will 
be able to address this question. 
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